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Abstract 



The cosmological constant (CC) problem is the biggest enigma of theoretical physics 
ever. In recent times, it has been rephrased as the dark energy (DE) problem in order to 
encompass a wider spectrum of possibilities. It is, in any case, a polyhedric puzzle with many 
faces, including the cosmic coincidence problem, i.e. why the density of matter pm is presently 
so close to the CC density p\. However, the oldest, toughest and most intriguing face of this 
polyhedron is the big CC problem, namely why the measured value of pA at present is so 
small as compared to any typical density scale existing in high energy physics, especially 
taking into account the many phase transitions that our Universe has undergone since the 
early times, including inflation. In this Letter, we propose to extend the field equations of 
General Relativity by including a class of invariant terms that automatically relax the value 
of the CC irrespective of the initial size of the vacuum energy in the early epochs. We show 
that, at late times, the Universe enters an eternal de Sitter stage mimicking a tiny positive 
cosmological constant. Thus, these models could be able to solve the big CC problem without 
fine-tuning and have also a bearing on the cosmic coincidence problem. Remarkably, they 
mimic the ACDM model to a large extent, but they still leave some characteristic imprints 
that should be testable in the next generation of experiments. 
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1 Introduction 



High Energy Physics is described by quantum field theory (QFT) and string theory. Unfortu- 
nately, these theoretical descriptions are plagued by large hierarchies of energy scales associated 
to the existence of many possible vacua. Such situation is at the root of the old and difficult CC 
problem [l], i.e., the formidable task of trying to understand the enormous ratio between the 
theoretical computation of the vacuum energy density and its observed value, p\ ~ 10"^'' GeV^, 
obtained from modern cosmological data [5. The extremal possibility occurs when the Planck 
mass Mp ~ 10^^ GeV is used as the fundamental scale; then the ratio Mp/ p\ becomes ~ 10^^^. 
One may think that physics at the Planck scale is not well under control and that this enormous 
ratio might be fictitious. However, consider the more modest scale v = 2 My// 9 — 250 GeV of 
the electroweak Standard Model (SM) of Particle Physics (the experimentally most successful 
QFT known to date), where Mw and g are the boson mass and SU{2) gauge coupling, 
respectively. In this case, that ratio reads \{V)\/p\ > 10^^, where {V) = —{l/8)M^v'^ < is 
the vacuum energy (i.e. the expectation value of the Higgs potential) and Mh > 114.4 GeV is the 
lower bound on the Higgs boson mass. Although one may envisage the possibility that there is a 
cancelation between the various theoretical contributions to the physical CC (including the bare 
value), this has never been considered a realistic option owing to the enormous fine-tuning that 
it entails (which, in addition, must be corrected order by order in perturbation theory). 

In this Letter, we discuss a dynamical mechanism that protects the Universe from any initial 
CC of arbitrary magnitude \pW ^ p^, which could emerge, for instance, from quantum zero- 
point energy (contributing roughly ~ for any mass m), phase transitions {p\ = {V)) or even 
vacuum energy at the end of inflation. We admit that p\ = pA{t) (with p\{ti) = /oj^, pA(to) = /^a) 
can actually be an effective quantity evolving with time. 

Phenomenological models with variable p\ have been considered in many places in the liter- 
ature and from different perspectives, see e.g. [3]. At the same time, models with variable CC 
with a closer relation to fundamental aspects of QFT have also been proposed [U [Sj [6l [7]. In all 
these cases, the effective quantity p\ = ph{t) still has an equation of state (EOS) = — pA and, 
in this sense, it can be called a CC term. 

The basic framework of our proposal is the generalized class of AXCDM models introduced in 
[8], in which there is a fixed or variable p\ term together with an additional "effective" component 
X (in general not related to a fundamental, e.g. scalar, field). This particular class of variable 
CC models is especially significant in that they could cure the cosmic coincidence problem [8] 
in full consistency with cosmological perturbations |9j. Here we present a generalization of these 
models that might even cure the old ("big") CC problem [Tj. Recently, in [TO] a model along these 
lines was introduced with a DE density pD and an inhomogeneous EOS po = ^ Pd — I3H~°^ 
which includes a term proportional to the negative power of the Hubble rate H . This additional 
term becomes sufficiently large to compensate an initial p\ when this is about to dominate the 
universe and forces it eventually into a final de Sitter era with a small CC. For recent related 
work on relaxation mechanisms, see e.g. [HI [121 [13] . In a different vein, the CC problem can also 
be addressed in quantum cosmology models of infiation, through the idea of multiuniverses [H] 
and the application of anthropic considerations . 

Let us recall that, historically, most of the models addressing the relaxation of the CC were 
based on dynamical adjustment mechanisms involving scalar field potentials [15]. In the present 
work, the relaxation mechanism that we propose is also dynamical, it does not require any fine- 
tuning and, as noted, it does not depend in general on scalar fields. To be more precise, the 
model we present here is a AXCDM relaxation model of the CC, which includes also matter 
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and radiation eras. We study the two possibilities p\ < and p\ > 0, with arbitrary value. 
For p\ < 0, our scenario avoids the big crunch at early times and allows the cosmos to evolve 
starting from a radiation regime with subsequent matter and de Sitter eras like the standard 
ACDM model. Finally, let us emphasize that our method to tackle the CC problem is formulated 
directly at the level of the (generalized) field equations, rather than from an effective action 
functional. In this sense, we follow the historical path of Einstein's derivation of the original field 
equations. At the moment, a version of our model with an action functional is not available, 
but its efficiency at the level of the field equations is truly remarkable, as we will show. In this 
sense, its phenomenological success may constitute a first significant step in the way of finding a 
solution of the difficult CC problem. 

The present Letter is organized as follows. In section [2] we present the basic setup of our 
model. In section [3] we present a toy model of the CC relaxation mechanism which helps to 
understand the basic idea behind our proposal, although it is still too simple to describe our 
Universe. Only in section [3] we present a first realistic version of the full relaxation mechanism 
and we perform a numerical analysis of it. In section [U we discuss in more detail some aspects 
and implications of our model. Finally, in the last section we draw our conclusions. 



2 The setup 

We start from the generalized Einstein field equations 

Rfiu - -j^QtJLvR = i^tll^ + T^u + QtJiy Pk,es) , (1) 

where is the energy-momentum tensor of ordinary matter - including the energy densities of 
radiation {pr) and baryons (ph). Furthermore, describes the X component {px), interacting 
with the effective CC term g^i, pA.eflf in such a way that the total density of the dark sector, 
Pd = PA,e{i + Px, is covariantly conserved (in accordance with the Bianchi identity). The effective 
CC density pA,e{i is given by 

PA,eS = Pa + Pinv ■ (2) 

Here, p\ is an arbitrarily large initial (and constant) cosmological term, and pi^v = Pmv{R, S,T) 
is some function of the general coordinate invariant terms 



R = R^,g^'' = 6H\l 



S = R^.uR'"' = 12H^ 



1 Y 3 

+ - 



2 V 4 



(3) 



which we have evaluated in the flat Friedmann-Robertson- Walker (FRW) metric in terms of the 
expansion rate H = a/a, and the deceleration parameter q = —da/a^ = —H/H^ — 1. We find it 
useful to write the structure of pinv in the form 

/'inv = J, (4) 

where /3 is a dimension 6 parameter and / = f{R,S,T) is a dimension 2 function of the afore- 
mentioned invariants. This form is particularly convenient since the function / must grow at high 
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energies and hence the vacuum energy is ultraviolet safe, i.e. in the early Universe p\rw and 
PA,eff ~^ Pa' where p\ is arbitrarily large but finite. 

The generalized field equations ([I]) fall into the metric-based category of extensions of General 
Relativity. However, at this point the following observation is in order. In the literature, the 
extensions of Einstein's field equations are usually of a restricted class, namely those that can be 
derived from effective gravitational actions of the form 



(fx 



M2 

-JLr + F{R,S,T) 

lOTT 



(5) 



This class of models may be called the 'T(R,S,T)-theories" as they are characterized by an ar- 
bitrary (albeit sufficiently differentiable) local function F of the invariants defined in Eqs. jS]), 
usually some polynomial of these invariants. Work along these lines has been put forward e.g. 
in [16]. The particular subclass of models in which the function F depends only on R, or "F(R)- 
theories", is well known and has been subject of major interest [T71 [TS] . 

However, as advertised in the introduction, in this work we formulate the relaxation mecha- 
nism directly in terms of the generalized field equations without investigating at the moment 
the eventual connection with an appropriate effective action. The reason is, basically, because we 
aim at maximal simplicity at the moment. To be sure, after many years of unsuccessful attempts, 
the CC problem has revealed itself as one of the most difficult problems (if not the most difficult 
one) of all theoretical physics; and we should not naively expect to shoot squarely at it and hope 
to hit the jackpot at the first trial, so to speak. In this sense, if we can find a way to solve, or at 
least to significantly improve, the problem directly at the level of the field equations, we might 
then find ourselves in a truly vantage point to subsequently attempt solving the CC problem at 
the level of some generalized form of the effective action of gravity. 

All in all, let us warn the reader that the connection between the two approaches (viz. the 
one based on the field equations and the functional one) is, if existent, non-trivial. In fact, we 
note that the presumed action behind the field equations ([I]) need not be of the local form ([5]), 
and in general we cannot exclude that it may involve some complicated contribution from non- 
local terms. These terms, however involved they might be, are nevertheless welcome and have 
been advocated in the recent literature as a possible solution to the dark energy problem from 
different perspectives [IS]. In the present work, we wish to put aside the discussion of these 
terms, and, for that matter, all issues related to the hypothetical action functional behind our 
field equations. Instead, we want to exclusively concentrate on the phenomenological possibilities 
that our framework can provide on relaxing the effective CC term ^ depending on the choice of 
the function /. In particular, if the initial p\ is a very large cosmological constant associated to 
a strong phase transition (e.g. some GUT phase transition triggering the process of inflation, the 
zero point energy of some field, or just the electroweak vacuum energy of the SM), the late time 
behavior of pA,efr can be sufficiently tamed (without fine-tuning) so as to be perfectly acceptable 
by the known cosmological data. 

The main aim of the present approach is thus of practical nature; if the CC problem can be 
efficiently tackled at the level of the field equations to start with (something that, to the best of 
our knowledge, has never been accomplished before), it should fit the bill as it can be already a 
crucial first step in the path to solve the CC problem - namely, before unleashing a more formal 
(and, predictably, even more difficult) theoretical assault to it at the effective action level. We 
leave this part of the investigation for future work, and we concentrate here on the potential 
phenomenological benefits of assuming a set of generalized field equations of the form ([T]). 
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3 Toy model 



Let us first illustrate the mechanism for a universe with only radiation and with / equal to just 
the Ricci scalar i?, so that Eq. jj]) becomes 

To start with, we consider the case p\ < (as in the SM case) and take (3 > 0. Let us also assume 
a spatially flat Universe. At late times, the de Sitter regime is realized and the deceleration 
approaches the value g — > — 1 while the Hubble rate becomes constant and very small, H ^ 
(of order of the current rate Hq). The vacuum energy density takes on the tiny observed value 
/OA,eff ^ Pa — Pai which is related to by the Priedmann equation: 3MpH^ = 87rp\. Since the 
initial value \pW is assumed to be much larger than p^, the final Hubble rate is approximately 
given by 

Hl^-^. (7) 

We observe that the large \pW is responsible for the small final value of H^,, provided the pa- 
rameter /3 has a suitable order of magnitude and without any need of fine-tuning because this 
relation does not include differences between large numbers (cf. [2U] for a discussion on fine-tuning 
issues). Moreover, this solution is stable. Indeed, the driving of to small values by the large 
and negative < becomes compensated by the positive second term in PA.efTj which grows 
as H decreases. On the other hand, a potential instability caused by a growing H would also be 
unharmful because it would make the second term decrease, so that p\ < Q would stabilize H 
again. 

At earlier times, H ^ H^, and the relaxation of p\ originates from the (1 — q) factor in the 
function /, namely the negative p\ drives dynamically the deceleration parameter q to larger 
values until q becomes very close to 1, which corresponds to radiation-like expansion. However, q 
cannot cross q = I from below since the (positive) second term in /3A,eff would dominate over p\ 
and stop the cosmic deceleration before q reaches 1. Summarizing, this simplified model keeps 
the enormous vacuum energy p\ under control at any time thanks to the relaxation mechanism 
implemented in the function / in Eq. Furthermore, it provides a reasonable expansion history 
with radiation- like expansion (g ~ 1, large) in the past and a stable de Sitter solution (g = — 1, 
H = < Hq tiny) in the future. The transition is smooth and happens when the Hubble rate 
is sufficiently small to ensure the CC relaxation as explained above. 



4 Full relaxation model 

The simple model discussed above is able to handle the large negative term p\ without abrupt 
changes in the expansion history. However, the model is unrealistic in that there is no matter era 
yet, because q will stay around the radiation domination value g ~ 1 until the de Sitter phase 
starts. Therefore, we have to make sure that the universe goes also through the matter epoch 
by completing the structure of pA,eS with a term proportional to (1/2 — q)~^, which would work 
like but with q = 1/2 as the stabilizing point for the next high H interval. For this purpose, 
we use the scalar invariants from Eq. A useful expression is to involve not only R but also 
5, as follows: 

R^-S = 24:H\2 - q) (1/2 - q) . (8) 
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Figure 1: General behavior of the deceleration parameter q in the relaxation CC model ([2]) (solid) 
and in the ACDM model (dashed) as a function of the cosmological redshift z. For any initial 
the universe goes through a radiation-dominated epoch (g = 1), a matter-dominated epoch 
[q = 1/2) and, eventually, into a final de Sitter phase {q = —1) with PA^efr — P\ ^ \Pa\- Concrete 
parameters in this plot are as in Fig. [2l 



Notice that this combination is proportional to (1/2 — q) and hence allows the relaxation of the 
vacuum energy in the matter era. Again, this expression alone would be unrealistic because 
it would enforce the Universe to linger in the matter era and would prevent the existence of a 
preceding radiation era. However, we can combine the three invariants R, S, T to finally form 
a much more realistic ansatz for / in Eq. i.e. in such a way that the radiation and matter 
epochs occur sequentially before the de Sitter universe is eventually reached in the infrared. 
Indeed, consider the expression 

/ = ^^ + yRT (9) 

= 4i?2ii_^K^+y.72i/6(i_^)(l + ^2)_ 
(1 - Q) 

We see that / is constructed such that the first term contains only two powers of H, which ensures 
that the cosmic evolution is reasonably close to that of the ACDM model during the matter and 
subsequent de Sitter stages. The second term contains T to provide a different scaling (~ H^) in 
terms of the expansion rate. As a result, the factor (1 — q) will dominate over (1/2 — q) for large 
values of H, i.e., during the radiation regime. 

The matter-radiation transition ("equality") happens when both terms in / are of the same 
magnitude and the corresponding time is fixed by the parameter y ~ H~g^, where Heg ~ lO^i/o is 
the Hubble rate at equality. The generic behavior of q for the relaxation model under consideration 
can be seen in Fig. [H The transitions between different epochs are not as smooth as in the ACDM 
cosmos, although this feature depends on the detailed form of the function /. For our illustrative 
purposes, the qualitative results given here should be sufficient to appreciate the virtues of this 
relaxation mechanism. Its main benefit is the complete insensitivity of the universe with respect 
to an initial cosmological constant p\ of arbitrary magnitude (which in the present example we 
have chosen to be negative). 
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Remarkably enough, our construction does not require to fine-tune any of the parameters 
of the model. Thus, the insurmountable problems associated to the traditional "cancelation" 
procedure can be solved automatically through this dynamical relaxation mechanism, which is 
triggered by pure gravitational physics (no scalar fields at all). The current value p\ of the effective 
vacuum energy and the corresponding low Hubble rate -ff* are fixed only by the magnitude of 
the parameter /3, which is the power of a mass scale M. Since \pW » p^, Eqs. ([4]) and ^ 
indicate that 

|/3| ^ m6 = \p\\ . f ~ \p\\Hl (10) 

where ^ Hq 10"^^ QgV. Remarkably, M can be of the order of a typical Particle Physics 
scale; if e.g. the initial vacuum energy is ~ M|,, then M < 100 MeV (i.e., of the order of the 
characteristic QCD scale Aqcd where the lowest phase transition occurs in the SM). The above 
relation (fTO]l can be rephrased in another suggestive way. Since g ~ — 1 in the eventual de Sitter 
regime - which starts approximately near our time -, we find that the current value of the CC 
(which is of the order of the asymptotic value p^eff) roughly given by the appropriate ratio of 
the two order parameters characterizing the most extreme phase transitions ever occurred in our 
Universe: 

PA - PA,efr - • 

So far, the discussion of the CC relaxation was based only on the form of pA,eS in Eq. jj]), 
and probably it can be implemented in various ways without losing its benefits (e.g. models with 
inhomogeneous EOS, modified gravity Lagrangian). In the following, we will discuss the concrete 
dynamics in a AXCDM-like setup [8], where the total energy density ptot = Pr + Pb + Pd includes 
the usual components of the known universe, pr and ph (i.e. radiation and baryons) as well as 
the extra contributions from the dark sector: po = Px + PA,eS- The conventional components 
are covariantly conserved leading to the usual scaling laws pr oc and pb oc a~'^. Since the 
dark sector does not interact with the conventional components, pD is conserved, too. Prom the 
Bianchi identity satisfied by the terms on the l.h.s. of Eq. ([T]) and the covariant conservation law 
of ordinary matter (V^T^™ = 0), the corresponding covariant conservation in the dark sector 
reads 

V^[T^^.+5M.PA,efr] =0. (12) 

Let us assume that X is a pressureless component. Computing the previous expression in the 
FLRW metric, it boils down to 

PA,cS + PX + 3Hpx = 0. (13) 

This equation shows that the two components of po are actually interacting. The fact that the 
EOS parameter of X is taken to be ux = (i.e. pressureless) is because X can then mimic (and 
can be referred to as) dark matter (DM). In this sense, the energy density of the dark sector, pD, 
can be thought of as being the sum of the DM and DE (interacting) energy densities, where the 
DE one is, in turn, the sum of the true cosmological constant p\ and the effective gravitational 
component pinv, i-e. Eq. ([2]). Thus, we have all the necessary ingredients to implement realistically 
our Universe within the relaxation model. 
The basic dynamical equations read: 

= =HiP^ = %pr+Pb + PX+PA,e^), (14) 

\^ / Pc Pc 

a 1 / 1 1 \ 

qH"^ = — = P"(^ + = (Pr + i^Pb + izpx - PA.cff , (15) 

^ Pc „ Pc \ ^ J 
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logjpCz+l) logio(z+l) 

Figure 2: Left: Normalized energy densities 0,n{z) = Pn{z) / ptot{z) for the CC, dark matter, 
radiation and baryons. The initial CC is p\ = — lO^'^p^!. Parameters at z = read: ^\ = 
0.73, 17° = 10"^, Vtl = 0.04, go = -0.6. The eventual de Sitter regime is ptot p\ - P^c- 
Right: Absolute energy densities Pn/ P^c of the CC, dark matter, radiation and baryons. Note 
that PA,cfr < and the plot shows |pA,eff|- 

where p° = ptot(io) = 3HqMp/{8it) is the current critical energy density. 

The various EOS parameters for n = r,b, X, A are Ur = 1/3 (accounting for photons and light 
neutrinos) and uJb = i^x = 0, wa = —1- Using these equations, in Fig. [T] we plot the numerical 
solution for q{z) and in Fig. [2] we solve for the normalized densities ^n{z) = pn{z)/ ptot{z)^ where 
we have assumed an initial cosmological constant p\ = — 10^°p|?. This should suffice to illustrate 
the great efficiency of this relaxation mechanism. Since y ~ 10^^ H^^, Eq. ^TOh yields /3 ~ 
(10~^ eV)®. Notice that in the presence of several phase transitions, the value of /3 = is fixed 
by the strongest one. We have seen above that M < Aqcd — 100 MeV for all transitions below 
the Planck scale. 

The following two points are in order concerning the role played by the X component. On the 
one hand, let us note that we have taken X as representing the full dark matter (DM) content of 
the Universe. This is a possibility, which we have chosen for definiteness in this presentation of 
the model, in part for simplicity and also because, then, the quantity po = px + PA,eS provides 
a kind of economical unification of the DM and DE parts into an overall, self- conserved, dark 
sector. However, this ansatz must be further elaborated. In particular, it must be confronted 
with the structure formation data from the analysis of cosmic perturbations in this model [21]. 
Another possibility would be, for instance, that the "conventional DM" is contained in what we 
have called the ordinary energy-momentum tensor in Eq. ([1]). In this alternative scenario, the 
total matter content is covariantly conserved, and X appears as a kind of additional entity in the 
DE sector, which would interact with the effective CC, and whose only purpose is to make the 
total DE covariantly conserved. Since we still have lox = 0, the entity X looks now more as a 
new form of DM that is integrated into the DE. 

On the other hand, once the dynamics of pA,eS is fixed, in this case through ^ and dS]), the 
evolution of the component X, in whatever of the two options discussed above, is completely 
determined by the local covariant conservation law (fT^ . This implies that X cannot be generally 
assimilated to a scalar field, because a dynamical scalar field with some particular potential has 
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Figure 3: Effective equation of state Weff and relative effective DE density JIde = PDE/ptot- 
Parameters are as in Fig. [2l 

its own dynamics. In this sense, X is to be viewed as an effective entity within the generalized 
field equations. In the particular option that we have analyzed in Fig.O it is supposed to mimic 
all effects associated to a real DM substratum. 

5 Discussion 

Let us now further describe the different stages of cosmic evolution in this framework. First, in 
the matter era, the total energy density is dominated by the DM component px rather than by 
the vacuum energy PA.eff- Due to g ~ 1/2 this can be understood easily by eliminating px from 
Eqs. (fTlTBll . 

- ^) = f (^^^ - ^^'^^) « " f + ■ ^''^ 

Therefore, this epoch behaves very similar to the ACDM matter era. Finally, the vacuum compo- 
nent becomes dominant at very late times and the Universe smoothly enters the eternal de Sitter 
regime with a very small positive pA,efr — Pa Pa- both eras, the relaxation model does not 
deviate much from the ACDM model, and the sign of the large initial vacuum density p\ is not 
relevant. 

Significant deviations from standard cosmology emerge in the radiation era, because there is 
no constraint that enforces PA,eS to be negligible. In fact, the exact behavior of px and PA.eff 
depends on initial conditions and the details of /. Nevertheless, in that epoch (for which g ~ 1 
and pb is negligible), subtraction of Eqs. (|14|15p leads to 

^ = ^(l-9) = ;^(2^x + 2pA,efrj. (17) 

Moreover, the second term in the function / in Eq. ((H) is dominant in the radiation era (by 
construction). Thus, while radiation dominates we have 

PA,eS = PA+y. 72i^6(l ^ g)(g2 + 1) = + ' ^' ^^^^ 
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and by eliminating the Ricci scalar R from the two previous equations, we obtain 

PA,eS=p\ + Jl ^ T- (19) 

With /3 ~ —p\Hq, y = {Heq)~'^ and the Hubble rate H(,q ~ lO^i^o at the radiation-matter 
transition we find that the variable 



7 = P^-^T7T^r^-P^-^^•(^)' (20) 



^0 P ^ ~Pa f 

y • 2AH^ ■ QHl ^ ^'^ ' 144 " V ^ 



becomes subdominant for very large Hubble rates H » H^q. Eq. (fT^ has two solutions for /OA,efr) 



/0± = g (^4pX -px± ^327 + (4pX + pxYj , (21) 

and the physical one has to be compatible with |/OA,eff| \pW at late times. Also, px is taken to 
be positive. 

For negative initial vacuum energy p\ < the following limits exist. At very high redshift, 
where px ^ ~^Pa^ effective vacuum energy acts like a true (and subdominant) constant, 

p+^^{ip\-px + \ipi + px\) =P\- (22) 

Thus, the universe evolves in a standard way, and the X component redshifts like non-interacting 
dust. At smaller redshift, in the range when ^ <C px ^ ~'^Pa holds, we find a (temporary) 
tracking regime: 

P+-1 -px- (4pX + Px))= -\px . (23) 

In this regime, PA^eff not only tracks the energy density of the X component, but also that of 
radiation p^. Indeed, in view of (f23]) . the conservation equation (fT3ll takes on the form 



Pi + mpi ~ (for both Pi = px,PA,es) (24) 
during this regime. Finally, for px ^ \Pa\ ^^e relaxation of the CC becomes obvious by 



p+^- ^4p\ -PX + m + px\ + + 0{{p\)-')j « \p\\. (25) 

Note that the analytical discussion is nicely supported by the numerical results shown in Fig. [2l 
This relaxation regime also exists for positive initial p\ > 0, when px ^ \Ai\^Pa^ 



U.J _ , _ ^ 327 , .,_2. 



p_ - ^4p\ -PX- m + px) - + Oiip\)-')j « p\. (26) 

Whereas for higher redshift, when px ^ V^bTi also find a tracking regime, p_ ~ —\px- 
However, differently from the p\ < case, this tracking behavior is persistent even for px ^ Pa- 
Consequently, at the end of reheating, the energy densities of radiation, dark matter and dark 
energy could be of the same order of magnitude. 

Because of the tracking relation px — — 4/9A,eff pr, we should care about bounds from 
nucleosynthesis. At that time {z < 10^), we have in our example = ^A + ^x ~ 3r2x/4 ~ 0.08 
versus Q.r ~ 0.9 (cf. Fig. [2]), and so the ratio fi/j/ilr ^ 10% is safe for standard Big Bang 
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nucleosynthesis (similarly as in [8]). On the other hand, this model offers also the possibility to 
solve the coincidence problem, in that px, Pr and /OA,eff are not very different during the tracking 
regime until the matter-radiation transition (cf. Fig. [2]). Note also that in comparison to ACDM, 
much more dark matter is allowed before the relaxation regime thereby weakening the constraints 
on dark matter properties. 

Finally, we discuss the effective EOS Weff, which is a useful tool for comparing interacting DE 
models with non-interacting ones. According to [22], Weff is given by the EOS of a self- conserved 
DE component /3de in an universe with the same Hubble rate H{z) and total energy density ptot 
as the relaxation model. Within this effective description, DM obeys the usual scaling law of 
dust px oc {z + Vf since the interaction with DE is absent. Thus, 

- 1, 1 + ^ 1 rfpDE(^) 
3 Pde(^) dz 

with /3de = Ptot — Px — Pr — Pb- Since Weff is more accessible for observations at low redshift, 
we magnify this range in the second plot of Fig. [3l In the relaxation regime, Weff follows mostly 
the EOS of the dominant component. Notice that Ode is well defined everywhere despite the 
behavior of c^eff- 



6 Conclusions 

In this Letter, we have addressed the old cosmological constant problem, i.e. the difficult problem 
of relaxing the value of the cosmological vacuum energy. The necessity of tempering this value 
ably and plausibly is absolutely crucial for a realistic cosmological evolution from the early times 
till today. Indeed, the vacuum energy of the early Universe is expected to be huge in Particle 
Physics standards, since the expansion history must drive through a series of phase transitions 
of diverse nature; in particular, it goes through a process of fast inflation (presumably associated 
to some Grand Unified Theory) and also through the spontaneous breaking of the electroweak 
symmetry. Finally, it undergoes the more modest chiral symmetry breaking transition, which 
occurs at the milder scale h-qcD = 0(0.1) GeV and is connected to the quark-gluon-hadron 
transition. Even if the latter would have been the only phase transition ever occurred, and the 
associated vacuum energy density would still persist, it would be a disaster for our Universe. The 
reason is that the value of the cosmological constant associated to that energy density would have 
accelerated the expansion history to the point of preventing the formation of any of the structures 
that we now see in our cosmos, as they would have been ripped off by the fast expansion rate 
during the first stages of formation. However, being the Standard Model (SM) of Particle Physics 
such an extraordinary successful theoretical and experimental framework, we must conclude that 
both of its important contributions to the vacuum energy (associated to the electroweak and 
strong interactions) must have been relaxed very fast after the corresponding phase transitions 
occurred; namely, sufficiently fast as to insure not only the possibility to form structures in the 
late Universe, but also to leave the nucleosynthesis process fully unscathed after the first minute 
of expansion. 

In this Letter, we have tackled a possible cure to this longstanding problem; we have pro- 
posed a dynamical relaxation mechanism of the vacuum energy that operates at the level of the 
generalized gravitational field equations. Our relaxation mechanism achieves this goal without 
fine-tuning. Apart from the ordinary baryonic matter, which in our model stays covariantly 
conserved, we assume that the dark sector is made out of an effective cosmological term pA,efr 
and another dynamical component X exchanging energy with it. The cosmological term of our 
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model is actually an effective one, in the sense that it is defined as the sum of an arbitrarily large 
cosmological constant, and a particular combination of curvature invariants, pinv, such that 
the sum behaves as an overall CC term p\^es{t) = p\ + P\m{i), but one that evolves with time. 
As we have said, the component X interacts with the variable PK,eS, but the total density of the 
dark sector, pn = Px + /OA,eff, is covariantly conserved. This construction fits into the class of 
the so-called AXCDM models existing in the literature [8j . Interestingly enough, in the present 
context, the X component can be interpreted as the dark matter (DM) and the contribution pinv 
can be viewed as a (dynamical) dark energy component that adds up to the traditional cosmo- 
logical constant term. However, such dynamical DE component has nothing to do with scalar 
fields as it is of purely gravitational origin. Therefore, the total energy density of the dark sector 
splits into the sum of the various components DM, CC and DE, namely pn = Px + p\ + Pinv, 
whereas the ordinary (baryonic) matter does not interact at all with the dark sector and remains 
safely conserved. 

The evolution of this AXCDM universe keeps the total DE sub-dominant during the radiation 
and matter epochs, and only at late times it leads to a tiny effective CC whose smallness is a 
direct consequence of the large initial vacuum energy, p\, rather than to a severe fine-tuning 
involving ugly cancelation of large terms. The resulting cosmos can transit from a fast early 
inflationary Universe, then drive through the standard radiation and matter dominated epochs 
and, eventually, ends up in an extremely slow de Sitter phase; in fact, a model of Universe very 
close to the one suggested by the modern cosmological data [2]. Finally, we obtained new insights 
into the coincidence problem, as we observed an interesting tracking behavior in the radiation 
era. A more detailed exposition of our approach discussing the universality of the CC relaxation 
and the corresponding analysis of the cosmological perturbations will be presented elsewhere [21]. 
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